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D/H in a new Lyman limit absorption system at z = 3.256 
towards PKS1937-1009 
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ABSTRACT 

We have identified a new Lyman Hmit absorption system towards PKS1937-1009, 
with logiV(Hi) = 18.25 ± 0.02 at z = 3.256 that is suitable for measuring D/H. 
We find a 68.3% confidence range for D/H of l.Qt°^'§^ x 10^^^ and a 95.4% range of 
1.6l[J j X 10"^ The metallicity of the cloud where D/H was measured is low, [Si/H] 
= —2.0 ±0.5. At these metallicities we expect that D/H will be close to the primordial 
value. Our D/H is lower than the D/H value predicted using the calculated from 
the cosmic background radiation measured by WMAP, 2.60;'Iq ]^7 ^ 10~^. Our result 
also exacerbates the scatter in D/H values around the mean primordial D/H. 



Key words: nuclear reactions, nucleosynthesis, abundances 
lines - cosmological parameters. 
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1 INTRODUCTION 

. There are interesting disagreements between the various 
' ways of calculating the baryon density of the universe (O;,). 

The most precise Q,b measurement is from the recent WMAP 
■ obser vations of the CMB power spectrum dSoergel et alJ 
'|20o3)- Another way to measure Q,^ is by measuring the rela- 
tive abundances of the light elements produced in big bang 
nucleosyntheis (BBN), since in standard BEN theory Q,b is 
the only free parameter determining these abundances. At- 
tempts have been made to measure the primordial *He, ^Li 
and deuterium (D) abundances. In principle, Q.b is most reli- 
ably determined by D, due to its apparently uncomplicated 
evolution with time {Epstein, Lattimor & Schramm IQtJ) 
and because it depends more sensitively on Q.^ than ''He and 
^Li. There is broad agreement between the D and CMB 
Ob, but the ''He and ^Li measureme nts predict a lower Q.b 
\Coc et alJl2004l: ICvburt et alJl2003l) . Either there is some- 
thing wrong with our interpretation of the *He and ^Li mea- 
surements, or the standard BBN theory may need to be 
modified (e.g. seelBarger et aljliool iDmitriev et al.ll20ol 
llchikawa fc Kawasakill2004^ . 

The mean of the primordial D values as measured in 
quasar absorption systems is consistent with the WMAP 
CMB Q.b. However, there is statistically significant scatter 
of the individual D/H measurements about the mean. The 
metallicity of the absorption clouds used to make the D/H 
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measurements is low, [M/H]< —1.5, so that no astration (D 
destructio n due to stel lar nucleosynthesis) is thought to have 
occurred iFieldjll996l) . In addition, the primordial D/H is 
thought to be isotropic and homogeneous, so the scatter is 
hard to explain. One explanation is that the systematic er- 
rors in measuring D/H have been underestimated by the 
authors. However, with only half a dozen quasar absorp- 
tion system D/H measurements other explanations, such as 
some early mechanism for astration or a non standard BBN, 
cannot be ruled out. More D/H measurements are needed 
to resolve this problem. It is exciting to see further results 
being added to this short list: in the last year two D/H mea- 
surements have been p ublished: one toward s another QSO 
absorption line system (iKirkman et al .120031) and another in 
a low metallicity gas cloud near our Galaxy iSembach et alJ 
l2004l) . 

Here we present an analysis of D/H in a newly discov- 
ered z = 3.256 absorption syste m towards PKS 1937- 1009 
(see Fig.[Tl Fig.HJ. Coincidentallv. iBurles fc Tvtled (ll998a^ 
measured D/H in a different absorber at z = 3.572 towards 
this same QSO. The new absorber is not connected with the 
z = 3.572 absorber, and we expect it to give an independent 
measurement of D/H. 



2 OBSERVATIONS 

PKS 1937- 1009 was first identified 
Zem = 3.787 bv lLanzetta et all il99ll) 



as a quasar with 
Our observations of 
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Figure 1. The two dimensional LRES spectrum, showing the 
position of the z = 3.572 LL at ~ 4170 A. The LL due to the 
z = 3.256 system we are interested in can be seen at ~ 3880 A. 

PKS1937-1009 were taken on Keck I with the HIRES and 
LRES instruments. The LRES observations were taken on 
9th August, 1996 and the HIRES observations were taken 
over the 31st May and 1st June, 1997. The HIRES spec- 
trum was taken using the C5 decker, giving a resolution of 
8.8 kms"^ FWHM. The LRES spectra were taken using the 
300/5000 grating, giving a resolution of ~ 7.4 A FWHM. 

The data reduction was carried out using standard 
IRAF routines. The raw CCD frames were flat fielded and 
ADU counts corrected to photon counts. Medians were 
formed of the 2D images to eliminate most of the cosmic 
rays. An observation of a standard star was used to trace the 
spectra in the 2D CCD image (for both LRES and HIRES 
data) and that trace was adopted in carrying out the extrac- 
tion from 2D to ID spectra. The stellar trace was applied to 
the quasar exposures, allowing it to translate in the spatial 
direction to account for shifts between the relative placing 
of the quasar and star along the spectrograph slit. Opti- 
mal extraction was done, using the weights derived from 
the standard star spectrum. Error arrays were propagated 
throughout the extraction procedure. Polynomials were fit- 
ted to ThAr lamp exposures to relate pixel number to wave- 
length. We checked for shifts between calibration exposures 
bracketing the quasar observations and time-interpolated 
where necessary. The wavelength scales were corrected to 
a heliocentric reference frame and all spectra were rebinned 
onto a linear wavelength scale. 

The LRES spectra were flux calibrated using the stan- 
dard star observations to recover the correct underlying 
spectral shape of the quasar continuum. No flux calibration 
was applied to the HIRES spectra. The final HIRES spec- 
trum covers wavelengths from 4100 — 6400 A with a signal 
to noise ratio (S/N) per pixel of 30 at 5000 A. The LRES 
spectrum covers the range 3700 — 7400 A, with a S/N per 
pixel of ~ 200 at 6000 A. 



3 ANALYSIS 

We fit Voigt profiles to the relevant absorption lines using 
the program VPFIT^. This uses a minimisation proce- 
dure to find the best fitting redshift, column density (cm~^) 
and b parameter (by definition, h — \/2a where a is the 
Gaussian velocity dispersion) for each absorption line. The 
diagonals of the covariance matrix give the error estimates 

^ Carswell et al., http//www. ast.cam.ac.uk/~rfc/vpfit. html 



on each of the fitted parameters. These errors assume that 
the parameter space around the best fitting solution is 
parabolic, which is not always the case. In particular, if there 
are degeneracies in the fitted parameters, the parameter 
space needs to be explicitly explored to find the correct er- 
ror ranges. For the majority of the transitions we fit in the 
z = 3.256 system, the errors given by the covariance matrix 
are sufficiently accurate. However, to find the errors on the 
best fitting value of D/H, we do not use the errors given by 
VPFIT, but instead explicitly calculate as a function of 
D/H, as described in section 

3.1 Absorption lines present in the z = 3.256 
system 

Cll (1334), Sin (1260) and Siiv (1393 and 1402) (Fig.l^l, Hi 
Lya and Ly/3 (Fig.|lJ transitions are present at z = 3.256 in 
the HIRES spectrum. There is a narrow line in the blue wing 
of the Hi Lya line which we believe is D I. The z = 3.256 
Lyman limit (LL) is present in the LRES spectrum (Fig.|21l. 
All higher order Hi transitions are at wavelengths shorter 
than the LL due to the higher redshift absorber at z — 3.572, 
and cannot be detected in our HIRES spectrum. 

The HIRES spectrum covers the position of Sim 
(1206.5) and Fell (1144), but they are blended with Hi ab- 
sorption in such a way that no useful constraints are possi- 
ble. 

Both the Si ll and C ll lines fall in the Lya forest, where 
there is the danger that an H I forest line may be misidenti- 
fied as a metal line. However, we are confident that we have 
correctly identified Si ll and C ll for two reasons. Firstly, their 
widths are much smaller than expected for the average Lya 
forest Hi line. Secondly, they fall within a few kms"'^ of 
the expected redshift based on z(Hl) and z(Silv). Unfortu- 
nately, both Si II and C ll lines are to some extent blended 
with Hi Lya absorption, which may be hiding weaker Sill 
and C II components. 

The Si IV lines are at wavelengths longer than the quasar 
Hi Lya emission and so are not blended with forest lines. 

The C II and Si ll appear to be at the same red- 
shift, z(Cii)= 3.256009 ±0.000006 and z(Siii)= 3.256001 ± 
0.000009. 

There are three Silv components, one main compo- 
nent at +3.4 ±0.4 kms~^(redward) relative to the Cll/Sill 
position, and two lower column density components at 
±21.4 ±0.7 kms"^ and ±43.4 ± 0.7 kms"^ 

The metal line parameters and their errors are given in 
Table □ 

3.2 Hi column density 

The z = 3.256 Lyman limit is optically thick, so there 
is a lower limit on the total log A'^(Hl) of ~ 17.6. There 
are also damping wings at the Lya line, which further 
constrain Af(Hl) (Fig. 01. There might be another expla- 
nation for these damping wings; perhaps the continuum 
dips over the Lya line, due to an emission line on either 
side of the Lya line . How ever, composite QSO spectra [e.g. 
IVanden Berk et alJ l|20^)] display no emission lines associ- 
ated with the quasar redshift at the relevant wavelengths 
to cause this apparent dip in the continuum. We believe our 
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Figure 2. The entire LRES spectrum (above) and an expanded region of the spectrum showing the LL 
due to the z = 3.256 absorber (below). The fit corresponding to model (2d) (see section f3.. Si is shown in the 
expanded spectrum. 



Table 1. The absorption line parameters and their errors, found using from the diagonals of the covariance 
matrix in VPFIT, for Cll, Sin and the three Si IV components. 



Ion 


log(Ar) (cm-2) 


z 


b ( kms-i) 


Cll 


13.600 ± 0.045 


3.256009 ± 0.000006 


6.9 ± 0.7 


Si 11 


12.875 ± 0.093 


3.256001 ± 0.000009 


5.0 ± 1.4 


Si IV 


12.820 ± 0.011 


3.256057 ± 0.000002 


6.5 ± 0.3 


Si IV 


12.138 ± 0.039 


3.256303 ± 0.000005 


1.3 ± 1.3 


Si IV 


11.970 ± 0.036 


3.256626 ± 0.000006 


2.8 ± 1.3 



two estimates of the continuum, described below, adequately 
account for any uncertainty in the shape of the continuum. 

3.3 Continuum placement 

To generate the continuum we fitted regions that were ap- 
parently free from absorption around each absorption fea- 
ture, with low order Chebyshev polynomials. The damping 
wings present in the H I Lya line provide the best constraint 
on the total A'^(Hl), but the derived 7V(Hl) value is sensi- 



tive to the placement and shape of the continuum above 
the damping wings. To take this into account, we fitted two 
different continua above the Lya line using 2nd and 3rd or- 
der Chebyshev polynomials. We also allow the continuum 
to vary above the Lya and Ly/3 lines by multiplying it by a 
parameter, initially set to 1, which is allowed to vary during 
the minimisation. We use two independent parameters 
to vary the continuum level above the Lya and Ly/3 lines. 
An example of the typical fitted values and errors for the 
two continuum parameters is given in Table |5| For all our 
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Figure 3. The Cn, Sin and Silv lines along with their best fitting models used to derive the parameters 
found in tabled The velocity zero point is the redshift of the Cll line, z = 3.256009. The histogram shows 
the data and the curves show the fitted model. The absorption lines away from the v = position in the top 
two panels are due to Lya forest lines. The three components in the bottom two panels are all due to Si IV. 
The normalised residuals ((data — model) /error) and their la ranges are shown below each spectrum, offset 
from the zero level and scaled for clarity. Horizontal lines at zero flux show the fitting regions. 



models, the fitted continua are within 5% of our initial guess 
for the continuum. 



3.4 H I velocity structure 

We try to determine the velocity structure in two ways: 
firstly, by fitting the available Hi transitions with as few 
components as needed to obtain an acceptable fit, and sec- 
ondly, by looking at the metal lines to indicate the possi- 
ble position of H I components. Ideally we would like to fit 
many higher order H I Lyman transitions to help determine 
the velocity structure, but only the Lya and Ly/3 lines are 
available. 

A single H I component cannot satisfactorily fit both the 
Lya and Ly/3 lines, so at least two components are present. 



We expect at least four hydrogen components, one at 
z(Cll/Sill) and three corresponding to the Silv component 
positions. However, we don't know what H I column density 
is associated with each one. We expect the C ll and Si ll com- 
ponents to have a greater proportion of their associated H 
gas in the form of H I, as their ionization potentials are close 
to that of Hi. Therefore, if there is no significant metallicity 
variation across the absorption complex, an Hi component 
near 2:(Cll/Sill) is likely to be the dominant component in 
the Lya and Ly/3 line. 



Below we describe the velocity models we explored. 
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Figure 4. The Lya and Ly/3 lines, along with a fit to data. The fit corresponds to model (2d), a two Hi 
component fit where the level of the continuum is allowed to vary above both the hya and Ly/3 lines. The 
continuum over the Lya line was generated using a 3rd order chebyshev polynomial. The velocity zero point 
is the redshift of the C II line. The histogram shows the data and the narrow curve shows the fitted model. 
The normalised residuals are shown beneath each graph. 



3.4.-1 Model (1): two Hi components, z(Hl) tied to z(Cll) 

We initially assemble a model using two Hi components, 
with the redshift of the bluer component tied to the C ll 
redshift. In this fit we include Dl components for both Hi 
components. We tie the redshifts of corresponding D I and 
H I components together, and fit a single D /H to both com- 
ponents (we discuss this assumption in section 13.81 . Since 
the redder D l component is heavily blended with the bluer 
Hi component, we tie its b parameter to the b parameter of 
the corresponding Hi component, assuming purely thermal 
broadening (this assumption has no effect on the final D/H 
value; see section [3.611 . We fit the summed total of the Hi 
column density of both components, instead of individually 
fitting the column densities of each component. We do this 
because while the column density in each component is not 



well constrained, the total log 7V(H l) is. This model does not 
fit the D I line in the Lya line adequately. We can obtain 
an acceptable fit if we include very weak absorption from 
some unknown contaminating metal line near the D I posi- 
tion. However, this results in D/H < 0.6 x 10^^ {la upper 
limit). This is smaller than even the lowes t estimates of the 
interstellar medium D/H, ~ 1.0 x lO"'^ jHebrard fc Mood 
l2003h . 

3.4-2 Model (2): two Hi components, z(H\) not tied to 
z (metals) 

It is possible that the main Hi component does not fall at 
exactly the same position as the C ll and Si ll lines. This 
could be due to several reasons: perhaps an error in the 
wavelength calibration of a few kms~^, or a gradient in the 
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Figure 5. The Lya line showing the damping wings and the best fitting model, model {2d). The histogram 
shows the data and the narrow curve shows the fitted model. The normalised residuals are shown below the 
spectrum. 



metallicity across the cloud. For this reason, we explore a 
model that does not tie the redshifts of the H I components 
to that of the metal lines. This model is identical to the 
model above, except it no longer ties the redshift of the 
bluer Hi component to the Cll line. This model gives an 
acceptable fit to the data. The position of the bluer hydrogen 
component is —3.5 ± 0.7 kms~^ relative to C ll. The redder 
Hi component is at +36 ± 4 kms~^. 

3.4-3 Model (3): four Hi components, z(Hl) tied to 
z (metals) 

We included four H l components, the redshift of each tied 
to a metal line component redshift (one corresponding to 
the Cll/Sill position and one for each Siiv component). We 
include Dl for all components where log A''(Hl)> 17.0. This 
does not give a satisfactory fit to the data: the narrow D I 
line in the Lya line in particular is very poorly fit. As in 
model (1), the fit could be improved by including very weak 
contaminating absorption near the D I position. Again, this 
gives a low D/H, < 0.5 x 10~^, much lower than the ISM 
D/H. 

Model (2) gives the only acceptable fit to the data. Mod- 
els (1) and (3) can give an acceptable fit, but we need to in- 
clude an unidentified metal line blended with D I, and there 
is no other evidence that this is present. If we do include such 
a fine, the measured D/H in this absorber is much smaller 
than the ISM D/H. 

We split model (2) into four sub-models, using (2a) a 
2nd order polynomial fixed continuum, (2b) a 3rd order 
polynomial fixed continuum, (2c) a 2nd order polynomial 
varying continuum and (2d) a 3rd order polynomial varying 
continuum. 



3.4-4 Two unresolved questions about the best-fitting 
model 

The best fitting model seems to fall below the data at —90 
kms~^ in the Ly/? line (Fig. [IJ. There could be several 
reasons for this. It is possible that the apparent discrepancy 



is due to noise, or there may be a weak, narrow Lya forest 
line near —75 kms~^. However, such a line would need 
to have a width much smaller than the median width of 
forest lines. A metal line could have a narrow width, but 
the higher redshift system at z = 3.572 does not have any 
metal lines which fall at this position. No other absorption 
systems were identified that could explain a metal line 
appearing at this position. Instead there may be a weak 
cosmic ray at this position that has not been completely 
removed. Two of the three 2D raw images combined to 
form the final spectrum show cosmic rays at the positions 
corresponding to —75 to —95 kms~^ in the Ly/3 line. We 
tentatively conclude that a noise spike is the cause of this 
discrepancy. We checked whether this discrepancy has a 
significant effect on the measured D/H by fitting the Lya, 
Ly/3 lines and LL excluding the region from —75 to —95 
kms~^. We found that the best fitting D/H value and its 
error are almost unchanged, as the tightest constraint on 
NCD l) comes from the Lya D I line. 

We would expect the Si ll and C ll redshifts to be consistent 
with the redshift of the Hi line associated with them. 
Instead we find that the Sill/Cll position differs from the 
Hi position by 3.5 ± 0.7 kms~^, or about 1.5 pixels. It 
is possible that there are problems with the wavelength 
calibration. This is unlikely, however, given that the Sill 
and Cll positions agree to within 1 kms~^, and that the 
typical wavelength accuracy of a polynomial fit to an arc 
exposure from HIRES is better than 0.5 kms~^. It is 
more likely that our velocity model is wrong in some way. 
The true velocity structure of the absorption complex will 
almost certainly be more complicated than our simple two 
component model. However, we note for all other velocity 
models that we considered, D/H was substantially lower 
(< 0.6 X 10~^) than the value that we quote for Model 2(b). 

Both of these questions may be resolved by much higher 
S /N observations of this absorption system. 
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Table 2. The best fitting parameters of the lines used to fit the Lya and Ly/3 lines and the LL of the z = 3.256 
absorber. The errors are those returned by VPFIT, calculated using the diagonals of the covariance matrix. 
The first two H I components in the table have associated D I lines which we included. We fitted the total 
column density of both these components - the column density entry next to the first H I component give 
the value of and error on the total column density of both Hi components. The bottom two rows give the 
best fitting continuum parameters (which are multiplied by the initial guess for the continuum level) for the 
Lyo and Ly/3 regions. 
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Hi 
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Hi 
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Hi 
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Hi 
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Hi 
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3.247846 
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0.000006 


21.5 ± 1.0 


Hi 


13.195 ± 0.073 


3.248529 


± 


0.000114 


63.0 ± 8.7 


Hi 


11.89 ± 0.23 


3.247192 


± 


0.000030 


7.7 ± 4.6 


Hi 


12.15 ± 0.11 


3.249996 


± 


0.000035 


14.8 ± 4.1 


Hi 


14.179 ± 0.031 


2.593103 


± 


0.000016 


32.6 ± 3.2 


Hi 


12.63 ± 0.38 


2.592337 


± 


0.000129 


20.3 ± 13.8 


Hi 


13.21 ± 0.11 


2.593948 


± 


0.000064 


27.9 ± 5.9 


Hi 


13.093 ± 0.094 


2.587331 


± 


0.000015 


12.1 ± 2.5 


Hi 


13.548 ± 0.038 


2.587413 


± 


0.000031 


53.4 ± 5.5 


Hi 


12.53 ± 0.10 


2.586028 


± 


0.000024 


11.2 ± 3.4 


cont{Lya) 


1.0035 ± 0.0027 










cont(Ly/3) 


0.983 ± 0.014 











3.5 Ionization and metallicity 

We have A'^(Sill) and Ai'(Cll) in a single gas cloud, 
and A^(Silv) in three other clouds. We use the program 
CLOUDY (^Ferland.l997il to generate an ionization model 
for the absorption system. When generating the model, we 
assume that the bluest Siiv component, the Sin and Cll 
component, and the bluest Hi component are all associated 
with the same gas, even though they span ~ 6 kms~^. Wo 
approximate the geometry of the cloud as a plane parallel 
slab illuminated by a distant point source and assume the 
proportion of metals to be solar. For calculating ionization 
fractions, the hydrogen particle density, uh and the inci- 
dent radiation flux on the cloud, Jv are degenerate. Instead 
of specifying either nn or J^, we use the parameter U: 

^ ^ 2 X io-.-/oi^/(10-"'^ergcm-^ s-^ sr"^) 
Uh /(lcm^'=') 

Here Jgi2 is the incident radiation at 912 A. We generate ion- 
ization models for a range of metallicities and U values, using 
a Haard t-Madau ionizing radiation continuum at a redshift 
of 3.25 faaardt fc Madaul [l99fi'). We have not explored the 
effect a different ionizing radiation continuum may have on 
the ionization fractions. A'^(Sill) and A'^(Silv) suggest that 
the system has a total metallicity [Si/H] ~ —2.0. Ai'(Cll) 
is lower than expected for this [Si/H] which might indicate 
that C is under abundant compared to Si (see Fig. |SJ . This 
is unlikely to be due to differential dust depletion, since C 



has a lower condensation temperature than Si. We would 
expect Si, not C, to be under ab undant if there is any deple- 
tion on to dust et al.llT996ri . However, an under abun- 
dance of C is consistent with metallicity measurements in 
metal poor stars and Hll regions, which find [C/O] ~ —0.5 
(Akcrman ct al. 2004). Since both Si and O are a-capture 
elements, it is reasonable to also expect [C/Si] ~ —0.5 when 
[Si/H] = —2 (Max Pettini, 2004, private communication). 

With so few transitions we can only make broad gen- 
eralisations about conditions in the cloud. We take a value 
of [Si/H] —2.0 with an error of ~ 0.5 dex, where this error 
is dominated by the uncertainty on Ai'(Hl) in the bluest Hi 
component. 

3.6 b parameter test 

We can predict the h parameter of the blue D I component 
from the b parameters of the Sill, Cll and Hi lines. This 
is done by measuring the thermal line broadening, fcthorm, 
and turbulent broadening, 6turb, by using the b parameters 
from all available ions. If we assume the absorbing cloud has 
a thermal Maxwell-Boltzmann distribution and any turbu- 
lence can be described by a Gaussian velocity distribution, 
then the b parameter for a particular ion will be given by 

^ion = ^thorm + ^turb ■ (1) 

Here fcthorm = where T is the temperature of the 
gas cloud, m is the mass of the absorbing ion and k is Boltz- 



8 




Figure 6. The ionization model generated using CLOUDY. The parameter U is defined in section [ij. 51 



mann's constant, fttu^b represents the Gaussian broadening 
due to small scale turbulence, and is the same for all ionic 
species. We plot 6^ against inverse ion mass in Fig.|7] Subject 
to the assumptions above, all the ions in the same cloud ve- 
locity space should lie on a straight line whose intercept gives 
feturb and slope gives the cloud temperature. C ll. Si ll and H I 
all have relatively low ionization potentials (IP); 24.4, 16.3 
and 13.6 eV, respectively. There is evidence from damped 
Lya absorbers that even 'intermediate' IP ions, such as Al ill 
(IP 2 8.4 eV), have a similar ve locity structure to lower IP 
ions (IWolfe fc Prochaskallioooll . Thus we assume that the 
Si II and C ll lines arise in the same gas as the H I and fit 
a least squares line of best fit to these three points. This is 
shown as the dashed line in Fig. |7| The D I point is plotted 
with its 1(7 error bars as given by VPFIT. We find a tem- 
perature of 1.35 ± 0.08 X lO^K and a turbulent broadening 
of 5.0 ± 0.8 kms~^. These values suggest a b parameter for 
the blue Dl component of 11.7 kms~^. This is consistent 
with our fitted &(Dl) of 11.6 ± 0.9 kms"^ for the blue Dl 
component. 

There is no way to measure the b parameter of the 
red Dl component, since it is heavily blended with the Hi 



blue component. We also cannot predict what b(Dl) will 
be using the method above, since we do not have any low 
ionization metal lines corresponding to the red Hi com- 
ponent. In our velocity models we explored the maximum 
and minimum 6(D l) values, corresponding to purely thermal 
[6(Dl)= I/V2 6(Hl)] and purely turbulent [6(Dl) = &(Hl)] 
broadening. The measured D/H is not affected by which 
broadening we assume. Thus we arbitrarily choose pure ther- 
mal broadening, assuming that the the conditions in the red 
component are similar to the blue component. 

3.7 Are we seeing deuterium? 

We expect we are seeing D, for three reasons. The line we 
have identified as D I occurs near the expected position for 
Dl, given the position of the main Hi absoprtion and the 
position of the Sill and Cll lines. The 6 parameter of the 
D I line is consistent with that expected for D I based on 
the fitted b parameters of the Hi component and the Sill 
and Cll lines. Finally, the line has a width substantially 
narrower than the median for forest lines, which means it 
is unlikely to be an Hi Lya forest line. A metal line from 
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Figure 7. The b parameter squared ( kms~^)^ versus the inverse 
ion mass (amu~^). The dashed Une is the least squares hne of best 
fit to the Si II, C II, and H I points. D I is shown for comparison. 
The error bars in each case are calculated from the la errors given 
by VPFIT. 



another absorption system could have a line width this nar- 
row, but we know of no other absorption systems towards 
this quasar which have metal lines that fall at the position 
of the D I line. These arguments all apply only to the bluer 
Dl component. We cannot independently measure the po- 
sition, h parameter or column density of the the redder D I 
component. However, we assume that D/H is the same in 
both components for the reasons given in section 



3.8 The validity of our assumption of constant 
D/H across different components 

We assume in our models that D /H is the same for all com- 
ponents showing D I. D is destroyed in star formation, so 
if the components have different metallicities, it is conceiv- 
able that they will also have different D/H values. Unfor- 
tunately, we cannot directly measure the metallicity of in- 
dividual components due to the Lya forest contamination 
around the Sin and Cll lines. However, we believe that 
even if the metallicity is significantly different in different 
components, D/H will still be very similar. Even in sys- 
tems where the metallicity has been found to vary consid- 
erab le between componen t s in the same absorptio n system 
(e.g. iKirkman et"ani2003l : iBurles fc TvtleI^ll998a^ . the dif- 
ference is ra r ely more than ~ 1.0 dex. T heoretical models 
jFieldjll99^ iPrantzos fc Ishimarull200 predict that D/H 
remains roughly constant with metallicity until metallicities 
approaching 1/5 solar are reached. Since [Si/H] ~ —2.0 for 
this entire absorption complex (see section \6.bl . we expect 
D/H in this absorber to be very close to the primordial value. 
Even if there are components with [Si/H] as high as ~ —1.0 



present, we don't expect their D/H to be significantly lower 
than that of any lower metallicity components. 



3.9 D/H in the z=3.256 absorber 

We have the values of A'^(H l) and A''(D l) and their errors for 
models (2a) - (2d), but we prefer to use an alternate method 
to find D/H in this absorber. Only two small regions con- 
strain A^(Dl), and the errors on A^( Dl) are relatively large. 
As we saw in lCrighton et aP (l2003li . when the errors given 
by VPFIT are large, it is often a sign that the parameter 
space around the best fitting solution is not symmetrical. In 
this case the space should be explored explicitly to find 
the correct error ranges. 

We have generated a list of values as a function of 
D/H. For each value of D/H, we fix D/H and fit the Lya, 
Ly/3 lines and the LL, varying all other parameters to min- 
imise . In this case we do not fit the the total Ai'(Hl) and 
A''(Dl), since we are not interested in the errors on individ- 
ual parameters, just the total value. In this fit we tie the 
redshifts of each Dl line to its corresponding Hi line. We 
also tie the h parameter of the red D I component to its cor- 
reponding Hi component, assuming thermal broadening, as 
discussed in section 

We calculate one of these graphs for an appropriate 
range of D/H values for each model, (2a) - (2d) (Fig. (HJ. 
The distribution of Ax^ = ~ Xmin for each graph is the 
same as that of a x^ distribution with a number of degrees 
of freedom equal to the number of fixed parameters (in this 
case, one - D/H). Here Xmin is the smallest value of x^ for 
a particular graph. For each point on each graph, we are 
fitting 990 data points with ~ 47 parameters, giving ~ 943 
degrees of freedom. The typical reduced x^ for our best fit- 
ting models is 1.2. 

D/H is more tightly constrained in the cases without 
a fioating continuum level, which is expected as there are 
fewer free parameters. The 3rd order continuum allows a 
better fit to the data than the 2nd, and is consistent with a 
larger D/H value. The predicted 68.3% range for D/H for the 



3rd order continuum models (2b and 2d) is 1.6^0 30 ^ 10 ^ ^ 
with a 95.4% range of l.et^ j x 10~^ The D/H ranges for the 
2nd order continuum models are lower, ~ 1.45 ±0.2 x 10~^. 
Model (2d), which takes the uncertainty in the continuum 
level [and any effect that may have on the uncertainty in the 
A'^(H l) value] into account, gives the best estimate of D /H in 
this absorber. It is worth noting that these error ranges are 
consistent with the ranges for D/H that would be derived 
by using the VPFIT best fitting values and errors given in 
Table. H 



4 COMPARISON WITH D/H MEASURED IN 
OTHER QSO ABSORBERS 

In this section we compare other D/H detections with 
our own detection, and the predicted D/H value from 
CMB and standard BBN. There are many upper limits 
on D/H in QSO absorbers that are much higher than the 
CMB predicted D/ H. Examples of such systems are to- 
wards Q 0014-H813 JCarswell et alJ Il994l : iRueers fc Hoear 
Il996allbl iBurles et al.l Il999^ . 00420-388 llcarswell et al " 
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Figure 8. D/H plots. The horizontal lines show the 68.3% and 95.4% confidence levels. The upper graphs are for models (2c) 

and (2d), where the continuum level is allowed to vary during the fit. The lower models are (2a) and (2b), where the continuum levels 
are fixed. The plots on the left use models with a 2nd order continuum fit over the Lya line, the plots on the right use a 3rd order 
continuum. 



19961 ), BR 1 202-0725 ('WamDier et al. ' '19961) , PG 1718- 
4801 JWebb et al. .1997; T vtlcr ct al. 1999; Kirkman et al 
' 200l|:ICrighton et alJl20oif and Q0130-4021 ijKirkman et al 



2000 ). We do not include these in our comparison, since they 
are clearly consistent with the predicted D/H value. 

We now describe six D/H detections in QSO absorbers, 
five of which we include in our comparison. We exclude one 
of the detections, towards Q0347-3819, for the reasons out- 
lined below in section 



the predicted b found using the Sill, Cll and Hi b pa- 
rameters. This is a grey LL system. A''(Hl) was measured 
usi ng the drop in flu x at the LL and initially estimated 
by iTvtler et al.l lll996l) as log NjUi) = 17.94 ± 0.05. This 
was revised inlBurles fc TvtleJ lll997l) to 17.86 ± 0.02 after 
IWampleJ 91996) and Songail a et al T Tl997ft pointed out that 
the continuum level of Tytler et al. was too high. The final 
D /H value from lBurles fc Tvtleil ^1998a^ is (3.3±0.3) x 10"^ 



4.1 z = 3.572 towards PKS1937-1009 

iBurles fc Tvtleil 1^213) have presented the most compre- 
hensive analysis of this system. Absorption attributed to D I 
is observed in two transitions. The measured 6(D l) matches 



4.2 z = 2.504 towards 1009+2956 

This system was analysed bv iBurles fc Tvtled ([19983). It 
has a grey LL and log A''(Hl) — 17.39 ± 0.06, measured us- 
ing the drop in fiux at the LL. D I is seen in only a single 
transition. The authors believe it is likely some contamina- 
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Figure 9. D/H measurements as a function of logAf(Hl) in the 
absorber wliere tlie measurement was made. The filled circle is 
the D/H value from this paper. The has hed region is the CMB + 
standard BBN prediction for D/H from lCoc et all bo04l . Refer- 
ences for the other D/H measurements are given in section Bl 



tion is present, so the Dl parameters are less certain than 
in other systems. As in our paper, the authors assume a 
constant D/H across the entire absorption complex. D/H 

= (4.ol:°-j) X 10-^ 



4.3 z = 2.536 towards HS0105+1619 

lO'Meara et al.l i200 J) analysed this system and found that 
6(Hl) is very small (~ 14 kms~^), meaning that each D I line 
is clearly separated from the nearby Hi absorption. Conse- 
quently 6(D l) and A''(D l) can be measured independently of 
the Hi absorption parameters, log iV(Hl) = 19.422 ± 0.009 
is determined by the damping wings in the Lya line, and D I 
is seen in five transitions. A relatively precise D/H can be 
derived: D/H = (2.54 ± 0.23) x 10-^ 

4.4 z = 2.076 towards QSO 2206-199 

iPettini fc BowenI J200 j) find a verv low fefHl) (~ 14 kms"^) 
and an extremely simple velocity structure (apparently a 
single Hi component). This allows a D/H measurement de- 
spite the relatively low S/N of their HST spectra. Two Ly- 
man transitions can be used to constrain A''(Dl). log Ai'(Hl) 
= 20.436 ± 0.008 is measured using the damping wings of 
the Hi Lya fine. D/H = (1.65 ± 0.35) x 10"^ 



4.5 z = 3.025 towards Q0347-3819 

This system was initially analysed by iD'Odorico et al.l 
l|200l|) . who found D/H = (2.24±0.67) x lO^It was subse- 
quently re-analysed bv lLevshakov et al.l ||20o3), who found 
D/H = (3.75 ± 0.25) x 10"^ log iV(Hi) = 20.626 ± 0.005 



is measured using the damping wings in the Lya line. Lev- 
shakov et al.'s result is different because they assume the 
main component in Hi has a similar z and 6 to that of 
H2, also observed in this absorber. The Dl lines are heavily 
blended with H I and neither &(H l) nor 6(D l) can be mea- 
sured directly. The D/H value derived is dependent on the 
assumed velocity structure, as shown by the different results 
of D'Odorico et al. and Levshako v et al. The Si iv and 4 tran- 
sitions for this absorber iu iProchaska fc Wolfe (1£99J show 
components at ~ —60 kms~^ (on the scale in Levshakov et 
al.'s Fig. 10), which may have associated Hi absorption that 
could affect the derived D/H. We suspect that the range of 
D/H values allowed in this absorber is significantly larger 
than that quoted bv lLevshakov et alJ Ill2002 ). If we were to 
include this result, we would choose a conservative Icr error 
range of (1.57 — 4.00) x 10~^. Since this range is at least 
a factor of 2 greater than the other D/H error ranges, we 
exclude this value from our comparison of D /H values below. 



4.6 z = 2.526 towards Q1243-)-3047 

iKirkman et al.l (|20o3) analysed this system. D I is seen in 
five transitions, but blended with Hi. A'^(Hl) is measured 
using the damping wings at Lya. This is complicated by 
the Lya lying close to the quasar H I Lya emission, where 
the continuum is steep and difficult to determine. The 
authors generate many continua based on 'anchor' points 
in the spectrum which they assume are free from forest 
absorption. They use these to find a robust estimate for 
log iV(Hl) = 19.73 ± 0.04. The fitted 6(Dl) is marginally 
consistent with the predicted value from 6(Hl) and b(0l). 
D/H = (2.42+°J^) X 10"^ 

All of these D/H measurements are in absorbers with 
metallicities of [Si/H] < —1, and so no significant D 
astration should have taken place. Thus we might expect all 
of these values, along with the D/H value from this paper, 
to be consistent with a single, primordial D/H. To test this, 
we fitted the five adopted D/H values above together with 
the value in this paper to a single D/H value. Where there 
are not symmetrical errors around the D/H values, we force 
them to be symmetrical by increasing the smaller error to 
match the larger error. We find D/H = (2.40 ± 0.3) x 10"^, 
which is consistent with the predicte d D/H from CM B and 
standard BBN, {2mto^^) x lO"'^ lICoc et aljl20o3) . The 
value of the total for the fit is 25.95 for five degrees 
of freedom. The chance of exceeding this value by 
chance is less than 0.01 %. Thus, taken at face value, the 
measurements are not consistent with a single D/H; but 
with so few measur ements this may just b e a result of small 
number statistics. iKirkman et alJ i2003l) argue that the 
inconsistency is probably due to underestimated systematic 
errors - indeed, if we artificially double all the errors of 
the D/H values, the total = 6.5 and the inconsistency 
disappears. However, we note that the statistical errors 
on each D/H result were arrived at only after careful 
consideration of systematic errors. This highlights the need 
for further D/H measurements to determine whether the 
scatter in D/H is real. 
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4.7 Correlation between D/H and iV(H l) 

IPettini fc BowerJ J200ll) and iKirkman et~ai] ll2003f) note 
that there appears to be a trend of decreasing D/H with 
increasing A'^(Hl). Our result does not follow this trend, and 
has a much lower D/H value than either of the two other 
low iV(Hl) systems towards PKS1937-1009 and 1009+2956 
fF i^. 1^. We are not sure of the reason for this difference. 
As lPettini fc BowenI i200 J) commented, it is interesting that 
the two systems with the highest D/H values are both grey 
LL systems, where the total Ai'(Hl) was measured using the 
drop in flux at the LL. In all the other D/H systems, includ- 
ing the one in this paper, the dominant constraint on the 
A'^(Hl) value is from the damping wings of the Lya or Ly/3 
lines. At first glance this might suggest a systematic error in 
measuring Af(Hl) could explain the difference between our 
D/H and the grey LL D/H values, although previous au- 
thors' analyses of systematics affecting A'^(Hl) in both grey 
LL and higher Ai'(Hl) systems have been very thorough. 



5 DISCUSSION 

We have identified a new Lyman limit absorption system 
towards PKS1937-1009, with logiV(Hl) = 18.25 ±0.02 at 
z — 3.256. The Lya and Ly/3 transitions are suitable for 
measuring D/H, and we find a 68.3% range for D/H of 
IStofo X 10~^ and a 95.4% range of 1.6t.ol x 10"^ 

The metallicity of the cloud where D/H was measured 
is low, [Si/H] — —2.0 ± 0.5. At this metallicity we expect 
that the D/H value will be primordial. 

This value disagrees at the 99.4% level with the pre- 
dicted D/H value using the i^i, calculated from the WMAP 
cosmi c background radiation measurements, 2.60lg j7 x 
10-5 jCoc et a]l2004) . and it exacerbates the scatter in D /H 
measurements around the mean D/H value. It is consistent 
with the D /H values measured in the ISM using FUSE, (e.g. 
iHebrarcffc Moos .2003,) . However, it is not consistent with 
the idea of a single primordial D/H value which stays con- 
stant until the metallicity reaches ~ 1/10 solar, then begins 
to drop off, reaching the ISM D/H at solar metallicity. Some 
early mechanism for D astration may be the cause for the 
scatter in D/H values iFields et alj|20()lh . although in this 
case we would expect to see a correlation between D/H and 

[Si/H], which is not observed. 

Other authors iPettini fc BowenI 1200 it iKirkman et alJ 

|2003^ note that there is a trend of decreasing D/H with 
increeising A'^(Hl) of the absorbers in which they are 
measured. The D/H value in this paper does not follow this 
trend, suggesting that the trend is not real and only due to 
the low number of D/H measurements. 
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